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Abstract 

The recent experimental data on Coulomb breakup of the nucleus '^^Ne are interpreted in terms 
of deformation. The measured large one-neutron removal cross-section indicates that the ground 
state of ^^Ne is either s- or p-halo. The data can be most easily interpreted as the spin of the 
ground state being 3/2~ coming from either the Nilsson level [330 1/2] or [321 3/2] depending on 
the neutron separation energy Sn- However, the possibility of 1/2+ coming from [200 1/2] is not 
excluded. It is suggested that if the large ambiguity in the measured value of Sn of '^^Ne, 0.29ibl.64 
MeV, can be reduced by an order of magnitude, say to be ±100 keV, one may get a clear picture 
of the spin-parity of the halo ground state. 

PACS numbers: 21.60.Ev, 21.10.Pc, 21.10.Jx, 27.30.+t 
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Recent experimental data obtained by using radioactive ion beams reveal that the neutron 
numbers such as N=8, 20 and 28 are no longer magic numbers in some nuclei towards the 
neutron drip line. Those neutron-rich nuclei can well be interpreted as being deformed 
and are often called nuclei in the island of inversion. One-neutron removal cross sections 
of a very neutron-rich nucleus ^^Ne due to the Coulomb breakup are reported in [Ij] and 
found to provide an evidence of the soft El excitation. The analysis of the cross section in 
Ref. [l| is based on one-particle wave-functions in the spherical Woods-Saxon potential. It 
was concluded that the measured large one-neutron removal cross sections were consistent 
with a P3/2 or si/2 neutron halo when a spectroscopic factor considerably smaller than 
unity was introduced. The heaviest halo nucleus so far experimentally established is ^^C 
with the ground-state spin of 1/2+ , namely an s-halo nucleus. While the nucleus ^^C may be 
interpreted as a spherical halo nucleus, the data on ^^Ne in |l| may be most easily interpreted 
in terms of the deformation or the one-particle motion in the deformed mean-field, since the 
spin-parity of the halo ground-state cannot be 7/2^ which is expected for the 21st neutron for 



spherical shape. We note that the analysis of t 
^jMgis and i|Ali2, is successfully performed in 



le spectroscopic data on light mirror nuclei, 
]^2ivx&i3 aiLKjL 13XT.112, lo Duv^v^coDiuiij pciivjiiiicu 111 [2j lu tcrms of the deformed mean field using 
Nilsson orbits occupied by the 13th nucleon. In the present Communication we attempt the 
interpretation of the data in [1] based on a deformed mean-field, using one-particle neutron 
wave-functions obtained by taking properly into account the weak binding. 

The measured low excitation energies 
are consistent with the picture that the Ne isotope with these neutron-numbers lies inside 
the island of inversion. Moreover, the large Coulomb breakup cross section reported in [l|, 
which clearly indicates the halo nature of the ground state of ^^Ne, suggests the efficient 
contribution by s or p component of the 21st neutron in the deformed mean field. These 
observations invite us to carry out the analysis of the data on '^^Ne in terms of deformation. 
It should be also noted that the measured magnetic moment of the ground state of f2Mg2i, 
which has the same neutron-number N=21, is reported in Ref. Sj and is consistent with 
the interpretation of /'^=3/2^, while both neighboring even-even nuclei, ^^Mg and ^^Mg, are 
well interpreted as being deformed. However, since the neutron separation energy of \ps/[g2i 
is 2.22 MeV, a neutron halo is not expected. 

Low-lying states of odd-A medium-heavy deformed nuclei are in a good approximation 
expressed by one (quasi)particle moving in the deformed field produced by the even-even 
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core 2!]. One-particle picture works much better in deformed nuclei than in spherical nuclei, 
because the major part of the long-range residual interaction in the spherical mean field can 
be included in the deformed mean field. In contrast, when the spherical shell model is applied 
to nuclei in the island of inversion, the resultant wave-functions are not easy to be predicted 
due to the complicated configuration mixing. Furthermore, since harmonic-oscillator wave- 
functions are always used in the traditional shell model calculation, the applicability of the 
calculations to the description of halo nuclei can be questioned. 

We start with our formulation for spherical shape. The neutron bound-state wave- 
function is an eigenfunction of the Woods-Saxon potential with an energy eigenvalue e < 

1 



$W : ijm) = ^R%{e,r)[Y, ® Xi^.n. (1) 



where 

R^^jj{e,r) oc arhe{ar), (2) 

with 

hi{-iz) = ji{z) + ine{z) (3) 

«^ = -^ (4) 

/"<ir|fli'|,.(e,r)p = 1 (5) 

We choose to express one-particle neutron wave-functions in the continuum using real energy 
variable e > 

I : ijm) = (e, r) [Y, ® X1/2U (6) 

^cos{6ij) kji{kr) - sin{5ij) kni{kr)) [Y^ (g) Xi/2]jm (7) 



V n'TTk 

where fj, expresses the reduced mass and 



while the normalization is expressed as 

drE!^[^{e,r)R^[^{e\r) = 5{e-e'). (9) 

The B(E1) value of the El excitation from a bound one-particle (jb) level to a continuum 
one-particle (jc) level keeping the even-even core in the ground state as a spectator is written 
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as 



l^e//ej ['--UcJb, 1, ^, ^,UjJ I y '^"'^^if, [£b,r)Kj^ {ec,r 

(10) 

First, it is found that the Coulomb breakup by the soft El excitation of halo neutrons 
occurs far outside the nucleus ^^Ne, noting that no one-particle resonance with the relevant 
angular momentum is present in the low-energy continuum. In Fig. 1 we plot the 2^3/2 
radial wave-function with the eigenenergy eb=—300 keV, the continuum si/2 radial wave- 
function with ec=+60 keV, and the product of those two wave-functions multiplied by radial 
variable. The value of Ec is chosen as an example, because for eb=—300 keV the quantity 
in Eq. fllOl) reaches the maximum around ec=+QO keV Q|. From Fig. 1 it is seen that the 
major part of the soft El matrix-element comes from the region of r=8-40 fm, far outside of 
the core nucleus ^°Ne. A similar result is obtained also for soft El excitations of the halo p 
neutron to continuum d levels, though the Coulomb breakup cross-section is much smaller 
than that of the excitation to s levels. Both the tail shape and the amplitude of the wave 
function of the halo neutron are very important for the cross sections of Coulomb breakup, 
while the wave functions inside ^^Ne may hardly play a role. 

In the phenomena where the structure of deformed wave-functions plays a role, the ex- 
pression in Eq. (ITOl) is replaced by the reduced matrix-element in Eq. (4-91) of Ref. [2I 
together with Eq. (lA-67) of Ref. [7|. Continuum one-particle wave-functions in a deformed 
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mean- field are calculated in terms of eigenphase [8|, |9| . For given Ec and ^l'" values there are 
a number of independent one-particle wave-functions, the number of which is equal to that 
of eigenphases, and all possible contributions must be calculated and summed up. However, 
when soft El excitations should be estimated which produce such large Coulomb breakup 
cross-sections as those in |l|, we may conveniently use the halo one-neutron wave-function 
taken from the s or p component of deformed Nilsson levels (an approximation in terms of 
"spectroscopic factor"), while one-particle wave-functions in the continuum are estimated 
for the spherical part of the Woods-Saxon potential. 

In Fig. 2 the Nilsson diagram is shown, of which the parameters are approxim ately 
adjusted to the (^°Ne+n) system. The calculation was done in the same manner as in jlO| . 
At (3=0 the 1/7/2 one-neutron resonance is found at 2.40 MeV with the width of 0.224 MeV, 
while neither the 2^3/2 nor 2pi/2 resonance defined by the eigenphase formalism j^, 9] is 
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obtained. However, the complicated non-linear behavior of the [330 1/2] resonant level for 
0.1 < j3 < 0.2 in the continuum (denoted by a dotted curve) indicates that the resonant-like 
component with Q'^=l/2^ coming from 2^3/2 is present around the energy region. Indeed 
the 2p3/2 resonance lying lower than the I/7/2 one is found if we use a slightly more attractive 
Woods-Saxon potential. For the parameters used in Fig. 2 the Nilsson level which is to be 
occupied by the 21st neutron is [330 1/2] for 0.22 < f3 < 0.30, [202 3/2] for 0.30 < (3 < 0.40, 
[321 3/2] for 0.40 < (3 < 0.59 and [200 1/2] for (3 > 0.59. Varying the parameters of the 
one-body potential within a reasonable range, any Nilsson levels other than those four levels 
are hardly obtained for the 21st neutron. The spin-parity of the lowest state is V=3/2~ for 
[330 1/2] due to the decoupling parameter lying between —2 (for P3/2) and —4 (for the /7/2), 
3/2+ for [202 3/2], 3/2" for [321 3/2] and 1/2+ for [200 1/2]. Among those four Nilsson 
orbitals the [202 3/2] level is excluded as a candidate for the configuration of the ground 
state of '^^Ne, because the smallest orbital-angular-momentum in the wave function of [202 
3/2] is £=2, which makes very little halo. Examining Fig. 2 we may also note that the 
presence of ^^Ne inside the neutron drip line is possibly realized by the deformation which 
is created by Jahn- Teller effect due to the near degeneracy of I/7/2, 2^3/2 and 2pi/2 shells in 
the continuum for spherical shape. 

In Figs. 3(a) and 3(b) the probabilities of appreciable components of the [330 1/2] and 
[321 3/2] levels calculated at (3=0.3 and 0.5, respectively, are shown, while the channels of 
Pi/2 (only in [330 1/2]), P3/2, /5/2, /7/2, ^9/2 and hii/2 are included in the calculation. The 
radius of the Woods-Saxon potential is fixed, while the depth is adjusted so as to obtain 



respective Nilsson levels as eigenstates of the deformed potential. As shown in Refs. [111. Il2|. 
the p components in = 1/2^ and 3/2~ Nilsson levels increase as the binding energies 
approach zero, though the probabilities at zero energies depend on Nilsson levels. This is 
in contrast to the fact that the probability of the s component in ^7^^=1/2+ Nilsson levels 
becomes always unity as the binding energy approaches zero. For example, at 6:^=— 300 keV 
the probability of the P3/2 component in the [330 1/2] for (3=0.3 and [321 3/2] for (3=0.5 
levels is 0.5225 and 0.2534, respectively. For reference, at 6:^=— 300 keV the probability of 
the Si/2 component in the [200 1/2] level for (3=0.5 is 0.67. 

Now, for example, the shape of the radial wave-function of the P3/2 component of [330 1/2] 
at (3=0.3 is not the same as that of the bound 2^3/2 level, since the latter is an eigenstate of a 
given spherical potential while the former is not. Moreover, two different spherical potentials 
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lead to different radial wave- functions of Si/2 at a given energy in the continuum. These 
differences may induce a non-negligible change in the resulting B(E1) values, even after the 
normalization of the 2^3/2 wave-function is adjusted to be the same as the probability of 
the P3/2 component in [330 1/2]. In Fig. 4 we show the squared radial integral on the r.h.s. 
of Eq. (ITU]) as a function of the continuum si/2 energy which is calculated using the 
following two kinds of bound-state wave-functions; (i) the 2^3/2 wave function with the 
energy eigenvalue —300 keV for a spherical potential and a normalization of 0.5225; (ii) the 
P3/2 component of [330 1/2] which has the energy of —300 keV in the deformed potential with 
/3=0.3. The appreciable difference in dB(El)/dE values of the two cases appears only for 
small Be values, and the integration of the two curves over e^. up till 2 MeV gives a difference 
of about 15 percent. Generally speaking, if we replace the P3/2 component of the [330 1/2] or 
[321 3/2] level by the 2^3/2 wave function with the same neutron-binding and normalization, 
the El transition 2^3/2 Si/2 has a larger B(El)-value than that of J93/2 Si/2- "Larger" 
or "smaller" depends on whether the spherical (2^3/2) eigenstate lies energetically higher or 
lower than the one-particle level in the Nilsson diagram for the relevant value oi (3 7^ 0. The 
larger the energy difference, the larger difference in B(E1) values appears which comes from 
the different radial shape of the two wave-functions. _ 

Combining Figs. 3(a) and 3(b) with Fig. 2 of Ref. pj], for a given Sn value we may 
find the relevant Nilsson level and the spin-parity of the ground state of the halo nucleus 
^^Ne, which are consistent with the observed Coulomb breakup cross-section. After taking 
into account both possible ambiguities in the parameters of the one-body potential and the 
approximation in terms of "spectroscopic factor" of s or p neutrons, we may conclude at 
least the following; The ground state has F=l/2+ coming from the Nilsson level [200 1/2], 
if Sn is appreciably larger than 500 keV. In this case the relevant deformation is expected 
to be very large, j3 > 0.6. If Sn is smaller than 500 keV, the ground state can be p-wave halo 
and F=3/2^. If Sn is smaller than 200 keV, the relevant Nilsson level may be [321 3/2]. 
Otherwise, it is [330 1/2]. 

To improve the accuracy of the measured Sn value much better than the available one 
3, S'„=0.29±1.64 MeV, can clarify the F value when it is combined with the data in 
ij. On the other hand, for example, the measurement of the magnetic moment of ^^Ne 
may not clearly pin down the spin-parity, since the estimated magnetic moment coming 
from the N=21st neutron in a deformed core is anyway negative and lies in the range of 
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—0.4 < nl hn < — 1-0 for possible states with J'^=3/2 coming from either [330 1/2] or [321 
3/2] Nilsson orbits and those with /'^=l/2+ from [200 1/2]. 

In the present analysis the possibility of excited states of ^°Ne after the Coulomb breakup 
is not included, however, the related Coulomb breakup cross-section is relatively small and 
lies within the ambiguities in the model calculation. The possible many-body pair-correlation 
is not included either. This is because, first of all, the halo neutron wave-function extends so 
much beyond the core nucleus that it couples weakly to the pairing field provided by the well- 
bound core nucleons. Secondly, in the ground state of odd-N nuclei the relevant neutron 
single-particle energy lies very close to the Fermi level. Then, irrespective of the nature 
of one-particle orbits, the occupation probability of the doubly-degenerate neutron level 
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15| . Therefore, 



obtained from solving the HFB equation is approximately equal to 0.5 
the contribution of halo neutrons to the Coulomb breakup which is estimated in the present 
article is expected to work first approximation. 

In conclusion, it is shown that the observed large Coulomb breakup cross section of ^^Ne 

n 

in Ref. p| is interpreted most easily and simply in terms of p-wave neutron-halo together 
with the deformed core ^°Ne. The measurement of Sn with an accuracy much better than 
the presently available one will clarify the spin-parity of the ground state of ^^Ne. 
Fruitful discussions with Prof. T. Nakamura are much appreciated. 
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Figure captions 

Figure 1 : The integrand of the radial integral on the r.h.s. of Eq. (ITOl) as a function of 
radial variable r. Namely, the product (solid curve) of radial variable, the bound 
2p3/2 radial wave-function Mpg/al^ — — 300keV, r) (dashed curve) in unit of fm~^/^, 
and the continuum si/2 radial wave-function R[^^^^{e = +60keV, r) (dotted curve) in 
unit of fm~^/^ MeV~^/^. The unit of the ordinate is for the product, while the figures 
written along the ordinate can be used also for the radial wave-functions, -^2^3/2 
R^sy2^ respective units. The radial integration of the solid curve gives the radial 
matrix-element of the El transition, 2^3/2 — > si/2- The depth of the Woods-Saxon 
potential is —46.75 MeV for which the eigenstate of the 2^3/2 neutron is obtained at 
—300 keV, while the diffuseness and the radius used are 0.67 fm and 3.946 fm (for 
A=30), respectively. 

Figure 2 : Neutron one-particle levels in Woods-Saxon potentials as a function of quadrupole 
deformation parameter (3. The potential depth is approximately adjusted so that the 
21st neutron of the prolately deformed nucleus ^^Ne can be a halo neutron. The 
depth, the diffuseness and the radius of the potential are —39 MeV, 0.67 fm and 3.946 
fm (for A=30), respectively. Positive-parity levels are plotted by solid curves, while 
asymptotic quantum numbers [N A Vt] are denoted for bound levels. See the text 
for details. 

Figure 3 : (a) Calculated probabilities of the major components of the [330 1/2] level with 
/3=0.3 as a function of energy eigenvalue e^; (b) Calculated probabilities of the major 
components of the [321 3/2] level with /5=0.5 as a function of energy eigenvalue e^,. 
The potential depth is adjusted to obtain respective Sq, values as energy eigenvalues 
of the deformed potential. The diffuseness and the radius of the potential are 0.67 fm 
and 3.946 fm (for A=30). 

Figure 4 : Estimate of | / r dr R^^^^^{-?,QQk(N ,r) Rf^^^{ec,r) p on the r.h.s. of Eq. as 
a function of Sc- The solid curve is obtained by using the ^3/2 wave-function taken 
from the ^3/2 component of the [330 1/2] level, which is bound at —300 keV for /?=0.3. 
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The probability of the P3/2 component is 0.5225. The dotted curve is calculated using 
the 2p3/2 wave-function with the energy eigenvalue —300 keV, which is normalized to 
0.5225. The depth of the Woods-Saxon potential of the former is —38.39 MeV, while 
that of the latter is -46.75 MeV. 
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(a) Probabilities of major components of tine [330 1/2] level 
R = 3.946 fm (A = 30) p = 0.3 
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(b) Probabilities of major components of the [321 3/2] level 
R = 3.946 fm (A = 30) p = 0.5 



"7/2 



P3/2 



'11/2 



5/2 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

10-9 -8 -7 -6 -5 -4 -3 -2 -1 

8^ (MeV) 



CD 



60 d 



50 ^ 



40 ^ 



30 ^ 



CM 



^ 20 



o 

V 



_J \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ L 

Kc|r|b)|^ for bound P3/2 with eb = - 300 keV 

and continuum S1/2 with £c 



P3/2 component of [330 1/2] for p = 0.3 is used 
2p3/2 wave function normalized to 0.5225 is used 



10 ^ 








—I — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — r- 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

8c (MeV) 



